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I.  Introduction 

This  report  summarizes  the  progress  of  our  studies  of  resonance  enhanced 
multiphoton  ionization  (REMPI)  of  molecules  and  molecular  fragments  and  of  single¬ 
photon  ionization  of  these  species  by  coherent  VUV  radiation.  The  focus  of  these  studies 
has  been  to  provide  both  a  robust  description  of  significant  spectral  features  of  interest 
to  applications  of  this  technique  such  as  ultrasensitive  detection  of  trace  and  transient 
species  and  needed  insight  into  the  underlying  dynamics  of  these  spectra.  A  summary  of 
our  progress  in  this  area  is  given  in  the  section  below. 

In  the  second  part  of  this  report  we  review  the  progress  of  our  effort  to  extend  our 
studies  of  photoelectron  spectra  from  the  energy  domain,  where  the  premium  is  on  high¬ 
ly  resolved  molecular  ion  spectra,  to  the  ultrafast  time  domain  where  the  emphasis  is  on 
exploiting  femtosecond  energy-and  angle-resolved  photoelectron  spectra  for  real-time  map¬ 
ping  of  wavepacket  dynamics  in  small  molecular  systems.  Such  time- resolved  photoelectron 
spectra  can  also  provide  unprecedented  insight  into  the  evolution  of  the  associated  elec¬ 
tronic  structure.  The  results  reported  here,  in  fact,  represent  a  significant  advance  towards 
exploiting  the  full  potential  of  femtosecond  energy-  and  angle-resolved  photoelectron  spec¬ 
troscopy  as  a  probe  of  wavepackets.  Such  time-,  energy-,  and  angle- resolved  photoelectron 
spectra  should  be  particularly  useful  for  studies  of  wavepackets  in  the  important  cases  of 
avoided  crossings  (nonadiabatic  transitions)  and  large-amplitude  motion.  The  key  to  these 
applications  is  a  robust  description  of  the  underlying  photoionization  amplitudes. 

II.  Resonance  Enhanced  multiphoton  Ionization  Spectra  of  Molecules  and 
Molecular  Fragments 

1.  Background  and  Objectives 
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Resonance  enhanced  multiphoton  ionization  (REMPI)  utilizes  laser  radiation  to  pre¬ 
pare  a  molecule  in  an  excited  state  via  absorption  of  one  or  more  photons  and  to  subse¬ 
quently  ionize  that  level  before  it  decays.  A  remarkable  feature  of  REMPI  is  that  the  very 
narrow  bandwidth  of  laser  radiation  makes  it  possible  to  select  a  specific  rotational  level  in 
the  initial  (ground)  state  and  to  prepare  the  excited  state  of  a  species  of  interest  in  a  single 
rotational  level.  Thus,  by  suitable  choice  of  the  excitation  step,  it  is  possible  to  selectively 
ionize  a  species  of  interest  which  occurs  in  very  minor  concentrations,  without  ionizing 
any  other  species  that  may  be  present.  This  conversion  of  optical  selectivity  into  chemical 
selectivity  makes  REMPI  one  of  the  most  powerful  tools  for  ultrasensitive  detection  of  trace 
and  transient  species.  On  the  more  fundamental  side,  coupled  with  high-resolution  pho¬ 
toelectron  detection,  REMPI  provides  ion  rotational  distributions  for  ionization  of  single 
rotational  levels  of  excited  electronic  states.  Such  state-resolved  spectra  can  clearly  provide 
significant  insight  into  the  underlying  dynamics  of  molecular  photoionization.  Other  appli¬ 
cations  of  REMPI  include  its  use  for  studies  of  state-selective  chemistry  and  for  exploring 
excited-state  chemistry  at  a  quantum-state-specific  level.  Finally,  although  REMPI  has 
the  distinct  advantage  that  a  single  rotational  of  an  excited  electronic  state  is  ionized  and 
that  it  can  be  achieved  with  photon  energies  less  than  the  ionization  potential,  studies  of 
photoionization  of  jet-cooled  molecules  in  their  ground  electronic  states  by  a  single  pho¬ 
ton  of  coherent  vacuum  ultraviolet  (VUV)  radiation  are  significant  in  their  own  right  and 
highly  complementary  to  those  of  REMPI. 

The  key  objective  of  our  effort  has  been  to  carry  out  quantitative  studies  of  REM¬ 
PI  of  molecules  and  molecular  fragments,  as  well  as  of  single-photon  ionization  of  these 
species  by  coherent  VUV  radiation,  in  order  to  provide  a  robust  description  of  significant 
spectral  features  of  interest  in  related  experiments  and  needed  insight  into  the  underly¬ 
ing  dynamics  of  these  spectra.  Joint  theoretical  and  experimental  studies  of  these  ion 
rotational  distributions  which  are  being  widely  studied  by  the  zero-kinetic-energy  (ZEKE) 
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technique  have  been  a  principal  focus  of  our  effort.^  This  technique,  which  is  based  on 
the  detection  of  photoelectrons  resulting  from  pulsed-field  ionization  of  very  high  Rydberg 
states  lying  just  below  an  ion  threshold,  makes  it  possible  to  obtain  cation  distributions 
with  sub-wavenumber  resolution.  The  unprecedented  resolution  of  this  ZEKE  technique 
has  opened  up  entirely  new  vistas  in  studies  of  photoionization  dynamics,  ion  spectroscopy, 
and  state-selected  ion-molecule  reactions.  Emerging  applications  built  on  the  ultra-high 
resolution  of  this  technique  include  its  use  for  accurate  determination  of  thermochemically 
important  ionization  energies,  for  characterization  of  ion  rovibrational  level  structure  of 
large  organic  molecules,  of  elemental  clusters,  and  of  weakly  bound  molecular  complexes, 
for  probing  reactive  fragments,  and  for  pump-probe  photoelectron  studies  of  wavepacket 
dynamics.  This  surge  of  experimental  activity  in  ultra-high  resolution  studies  of  molecular 
photoelectron  spectra  continues  to  raise  new  theoretical  challenges  and  has  provided  the 
stimulus  for  several  of  our  collaborations  with  experimental  groups  in  North  America  and 
Europe. 

2.  Highlights  of  Accomplishments 

We  will  now  use  a  few  examples  to  highlight  the  progress  of  our  studies  of  ion  rotational 
distributions  from  REMPI  of  excited  electronic  states  of  molecules  and  from  single-photon 
ionization  of  ground  state  jet-cooled  molecules  by  coherent  VUV  and  extreme  ultraviolet 
(EUV)  radiation.  Although  we  will  not  discuss  the  numerical  and  computational  details 
of  the  procedures  used  to  obtain  the  molecular  photoelectron  orbitals  needed  in  these 
studies,  it  is  essential  to  recognize  that  the  quantitative  determination  of  the  ion  rotational 
distributions  of  interest  here  requires  the  use  of  molecular  photoelectron  orbitals  which 
correctly  incorporate  the  angular  momentum  coupling  present  in  these  wave  functions.  In 
contrast  to  atomic  photoelectron  orbitals,  molecular  photoelectron  orbitals  are  not  angular 
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momentum  eigenfunctions  but  contain  admixtures  of  angular  momenta.  This  coupling 
of  angular  momenta  is  brought  about  by  the  torques  associated  with  the  nonspherical 
potentials  of  the  molecular  ions  in  which  the  photoelectron  moves.  This  simply  reflects  the 
fact  that  as  the  photoelectron  collides  with  the  ion  core,  its  angular  momentum,  as  well  as 
that  of  the  molecular  ion,  changes.  Such  angular-momentum  changing  collisions  between 
the  photoelectron  and  the  ion  clearly  play  a  crucial  role  in  determining  ion  rotational 
distributions. 

In  our  studies  we  take  the  photoelectron  orbitals  to  be  solutions  of  a  one- 
electron  Schrodinger  equation  containing  the  Hartree-Fock  potential  of  the  molecular  ion. 


(-^  +^«on(r,i2)  -  y)V’*(r)  =  0,  (1) 

where  fc^/2  is  the  photoelectron  kinetic  energy.  We  obtain  these  orbitals  numerically 
using  an  iterative  procedure,  based  on  the  Schwinger  variational  principle,  to  solve  the 
integral  (Lippmann-Schwinger)  equation  associated  with  eq.  (1).  A  signiflcant  feature  of 
this  procedure  is  that  it  avoids  the  integration  of  the  coupled  integro-differential  equations 
which  normally  arise  in  the  numerical  solution  of  Schrodinger  equations  with  nonspherical 
potentials,  V(r,  6,4>).  The  computationally  intensive  steps  of  this  procedure,  which  is 
numerically  very  stable  even  at  the  low  photoelectron  energies  (ss  0.05  eV)  often  of  interest, 
are  single-center  expansions  of  a  large  number  of  multicenter  basis  functions  and  operators 
and  subsequent  use  of  these  expansions  in  numerical  quadratures  of  radial  integrals.^ 


(a)  Rotationally  Resolved  REMPI  Photoelectron  Spectra  for  H2  O 

We  have  completed  combined  theoretical  and  experimental  studies  of  rotationally 
resolved,  REMPI  photoelectron  spectra  for  low  rotational  levels  of  H2O  at  photelectron 
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energies  from  about  1  eV  to  2.5  eV  above  threshold.  These  were  the  first  non-ZEKE  rota- 
tionally  resolved  photoelectron  spectra  for  single  rotational  levels  of  a  polyatomic  molecule. 
While  ZEKE  photoelectron  spectroscopy  provides  sub- wavenumber  resolution  in  ion  rota¬ 
tional  distributions,  the  technique  is  restricted  to  the  threshold  region  and  does  not  permit 
studies  at  photoelectron  energies  away  from  threshold.  Studies  at  higher  photoelectron 
energies  can  display  dynamical  behavior  which  may  not  be  apparent  from  near-threshold 
studies.  This  is  particularly  true  of  features  that  may  arise  from  Cooper  minima. 

Spectra  were  measured  for  (2  -f  1')  REMPI  of  specific  rotational  levels  of  the  C  (0, 
0,  0)  state  of  H2O  using  a  magnetic  bottle  photoelectron  spectrometer.^  This  spectrometer 
provided  a  kinetic  energy  resolution  better  than  4  meV.  Comparison  of  the  measured  and 
calculated  ion  distributions  is  very  encouraging  and  reveals  extremely  nonatomiclike  pho¬ 
toionization  dynamics  and  the  influence  of  Cooper  minima  in  the  photoelectron  channels. 
We  illustrate  these  results  with  spectra  for  the  4i3  rotational  level  of  the  C  (3pai) 
state.  Fig.  1  shows  our  calculated  rotationally  resolved  photoelectron  spectra  along  with 
the  measured  (2  -|-  1')  REMPI  spectra  for  photoionization  of  this  level.  The  calculated 
spectrum  is  convoluted  with  a  Gaussian  detection  function  with  a  FWHM  of  42  cm~^. 
These  spectra  show  only  type  h  transitions  [AKa  =  even  (odd),  AKc  =  even  (odd)],  with 
the  most  intense  transition  being  4i3  <—  4i3.  On  the  basis  of  the  selection  rules  governing 
this  transition,®  i.e.,  AKa  +  ^  =  odd  and  AKa  +  AKc  =  even,  this  transition  must  arise 
from  odd  partial  waves  of  the  photoelectron  matrix  element.  A  single-center  expansion  of 
the  3pai  orbital  shows  that  it  has  14.5  %  s,  83.6  %  p,  and  1.9  %  d  character.  In  view 
of  the  84  %  p  (^  =  1)  character  of  the  Spai  orbital,  this  intense  4i3  <—  4i3  peak  in  the 
photoelectron  spectra  is  highly  molecular  in  origin  and  arises  from  strong  ^-mixing  in  the 
molecular  electronic  continuum.  Similar  behavior  is  also  clearly  seen  in  the  mixed  peak 
for  the  533  and  615  ionic  levels.  While  the  other  AKa  =  even  and  AKc  =  even  {£  =  odd) 
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transitions  such  as  the  831,  Sis,  and  431  levels  are  also  dominant,  they  are  mixed  with  other 
rotational  peaks  arising  from  even  £  components  of  the  photoelectron  matrix  element. 

To  provide  some  insight  into  the  dynamical  aspects  of  the  angular  momentum  transfer 
of  the  photoelectron  upon  ionization  in  this  system,  we  show  the  photoelectron  matrix 
elements  for  kinetic  energies  up  to  10  eV  in  Fig.  2.  The  behavior  of  these  photoionization 
matrix  elements  is  very  nonatomiclike.  Surprisingly,  multiple  Cooper  minima  are  also  seen 
in  every  continuum  channel  around  2  eV.  These  minima  are  present  in  the  d  and  f  waves 
of  the  A:ai  channel,  in  the  p,  d,  f,  and  g  waves  of  the  khi  channel,  and  in  the  f  and  g  waves 
of  the  kb2  channel.  These  Cooper  minima  strongly  deplete  the  contributions  of  the  even 
and  odd  waves  over  the  photoelectron  energy  range  of  interest.  In  the  spectra  of  Fig.  1, 
the  contributions  of  the  odd  (especially  p  and  f)  waves  are  much  stronger  than  those  of 
the  even  waves.  The  partial  wave  mixing  near  threshold  is  particularly  strong. 

(b)  REMPI  Photoelectron  Spectroscopy  of  Vanadium  Dimer 

The  elucidation  of  the  electronic  and  geometric  structure  of  transition  metal  clusters 
remains  a  challenging  problem.  Although  a  wealth  of  spectroscopic  information  has  been 
obtained  on  neutral  transition  metal  dimers,^  very  httle  is  known  about  trimeric  and  higher 
clusters.  Moreover,  with  regard  to  transition  metal  cluster  cations,  there  is  a  severe  lack 
of  gas-phase  spectroscopic  data,  even  for  diatomics. 

ZEKE  photoelectron  spectroscopy  is  a  powerful  technique  for  addressing  these  prob¬ 
lems.  Recently  Yang  et  al.®  exploited  this  technique  to  study  the  electronic  structure  of 
the  vanadium  dimer  cation.  Rotationally  resolved  ZEKE  spectra  were  obtained  by  (1  -|- 
1')  REMPI  of  the  A  ^IIu  excited  state  of  V2.  Striking  differences  were  observed  in  the 
ion  rotational  spectra  recorded  through  the  ^IIiu  spin-orbit  component  and  those  recorded 
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through  the  ^n2u  component,  where  a  wide  range  of  values  are  accessed,  as  opposed  to 
only  two  in  the  case  of  the  ^IIiu  component. 

To  provide  a  quantitative  description  of  the  underlying  dynamics  of  these  features, 
we  have  calculated  the  ZEKE  photoelectron  spectra  of  the  X  ground  state  of 
for  (1  +  1')  REMPI  of  the  A  state.®  Fig.  3  shows  calculated  and  measured  ZEKE 
photoelectron  spectra  for  the  A  state  of  V2  with  the  excitation  laser  tuned  to  the  R(7) 
line  in  the  A  ^  Hi u  <-  X  band.  Agreement  between  these  measured  and  calculated 

spectra  is  quite  encouraging.  Due  to  selection  rules,  the  photoelectron  spectra  of  Fig.  3 
for  the  A  ^IIiu  state  show  transitions  only  to  the  spin-orbit  component  of  the  ion. 

Furthermore,  since  on  the  basis  of  parity  selection  rules,®  only  the  e  parity  component  of 
the  J'  =  8  rotational  level  is  accessed  in  this  R(7)  excitation,  the  number  of  J'*'  rotational 
levels  for  the  X  state  of  is  significantly  reduced.  The  agreement  between  the 
measured  and  calculated  spectra  is  reasonable  except  for  the  6.5(f)  and  9.5(e)  peaks, 
where  our  calculated  intensities  are  somewhat  stronger  than  the  measured  values. 

Fig.  4  shows  the  calculated  and  measured  ZEKE  photoelectron  spectra  for  (1  -|-  1') 

REMPI  of  the  A  ^Iliu  state  of  V2  with  the  excitation  laser  tuned  to  the  R(6)  line  in 

the  A  ^n2u  <-  X  band.  For  this  A  ^n2u  spin-orbit  component,  both  the  X 

and  X  components  of  the  cation  are  accessible.  Agreement  with  the  measured  and 

2  ^ 

calculated  spectra  is  encouraging  except  for  the  peak  for  the  =  8.5(e)  and  6.5(f) 
levels  which  shows  a  weaker  intensity  in  the  calculated  spectra.  Since  the  X  spin- 

orbit  component  of  the  ground  state  contains  both  e  and  /  parities  for  each  rotational 
level,  which  are  assumed  to  have  equal  populations,  both  the  e  and  /  parities  of  the  J'  =  7 
rotational  level  of  the  A  ^n2u  component  are  excited.  A  wider  range  of  J"*"  rotational  levels 
would  hence  be  expected  to  arise  here  than  for  excitation  via  the  A  ^IIiu  component. 
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(c)  ZEKE  Photoelectron  Spectra  of  Sodium- Water  and  Sodium- Ammonia 
Clusters 

The  energetics,  structure,  dynamics,  and  spectroscopy  of  clusters  consisting  of  neutral 
or  ionic  metal  atoms  and  polar  solvent  molecules  axe  of  significant  interest  in  studies  of 
chemical  reactions,  electron  solvation,  corrosion,  and  a  wide  range  of  chemical  and  biolog¬ 
ical  processes.  Such  clusters,  composed  of  a  metal  complexed  with  a  variety  of  solvents, 
particularly  those  containing  alkali  metals,  are  interesting  targets  of  study  with  ZEKE-PFI 
techniques.  ZEKE-PFI  spectra  can  be  very  useful  in  elucidating  the  molecular  structure 
and  intermolecular  force  fields  of  these  small  clusters.  A  combination  of  experimental 
studies  and  ab  initio  calculations  can  be  very  effective  in  unraveling  these  ZEKE  spectra 
and,  hence,  in  obtaining  structural  data  from  these  spectra. 

We  have  carried  out  such  a  joint  study  of  the  ZEKE  photoelectron  spectra  of  small 
Na(H20)n  and  Na(NH3)n  clusters  with  Geoffrey  Blake  of  our  Department.  Experimentally, 
it  is  straightforward  to  produce  these  clusters  in  sufficient  quantities  to  detect  with  ZEKE- 
PFI,  and  the  wavelengths  required  for  their  spectroscopic  investigation  are  easily  generated 
with  dye  lasers.  We  have  concentrated  on  measurements  of  the  ZEKE  rotational  band 
contours  in  the  sodium- water  and  sodium-ammonia  complexes  and  on  calculations  of  these 
rotational  spectra.  While  agreement  between  the  measured  and  calculated  spectra  is  quite 
encouraging,  the  more  significant  implication  of  these  results  is  that  they  demonstrate  how 
essential  calculated  spectra  axe  in  accounting  for  the  main  features  of  the  spectra  and  their 
potential  for  elucidating  the  molecular  structure  of  the  neutral  and  ionic  clusters.^ 

Fig.  5  shows  our  measured  and  calculated  ZEKE  photoelectron  spectra  for  ionization 
of  the  3sai  orbital  of  the  X  state  of  NaH2  0  by  single-photon  coherent  UV  radiation.^ 
The  calculated  spectra  were  obtained  for  a  photoelectron  energy  of  50  meV  and  were 
convoluted  with  a  Gaussian  detector  function  with  a  full- width  at  half-maximum  (FWHM) 
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of  2.2  cm“^  The  rotational  constants  used  were  derived  from  the  calculated  geometry. 
These  are  A  =  14.157,  B  =  0.280,  and  C  =  0.275  cm  ^  for  NaH20  and  A"*"  =  14.15<',  B"*" 
=  0.309,  and  C'*'  =  0.296  cm“^  for  NaH20‘'‘.  Comparison  of  the  measured  ZEKE  spectra 
with  spectra  calculated  over  a  wide  temperature  range  suggests  a  rotational  temperature 
of  about  100  K  in  these  experiments.  The  spectra  in  Fig.  5  were  calculated  at  this 
temperature. 

Agreement  between  the  meastired  and  calculated  spectra  is  encouraging.  The  domi¬ 
nant  peak  in  these  spectra  belongs  to  P  branch  transitions  {AN  =  —  !)■  The  intensity  of 
this  peak  is  due  to  a  band  head.  Analysis  of  the  calculated  spectra  shows  that  lAA'i  =  0-4 
transitions  contribute  13%,  62%,  15%,  7%,  and  3%  of  the  total  cross  section,  respective¬ 
ly.  The  dominance  of  lAA”!  =  1  transitions  in  the  ZEKE  spectra  for  photoionization  of 
this  3sai  orbital  with  its  91%  s  character  is  surprising  and  very  nonatomiclike  since,  on 
the  basis  of  angular  momentum  considerations,  i.e.,  \AN\  <  ^  -|-  1,  the  s(^  =  0)  wave 
of  the  photoelectron  drives  \AN\  =  1  transitions.  This  s  -4  s  behavior  is  highly  unusu¬ 
al  for  photoionization  of  an  essentially  3s  orbital  of  sodium  where  atomiclike  propensity 
rules  would  suggest  a  strong  p-wave  photoelectron  matrix  element.  In  fact,  our  calcu¬ 
lated  photoelectron  matrix  elements  show  an  unusually  strong  s-wave  component.  Even 
partial  wave  components  of  the  photoionization  matrix  element  account  for  about  70% 
of  the  total  cross  section.  Our  calculated  rotationally  resolved  photoelectron  spectra  and 
the  underlying  dynamical  coefficients  account  quite  well  for  the  principal  features  of  these 
spectra.^ 

To  explore  the  utility  of  studies  of  this  kind  in  elucidating  the  molecular  structure  of 
these  neutral  and  charged  clusters,  it  is  instructive  to  see  how  changes  in  the  geometry 
influence  the  contour  of  the  calculated  ZEKE  spectra.  With  such  small  B  and  C  constants 
(the  A  constant  does  not  play  a  role  in  determining  the  “end-over-end”  rotational  spacing), 
a  slight  change  in  B  and  C  can  lead  to  dramatic  change  in  the  spectral  profile  since,  at 
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a  temperature  of  100  K,  more  than  300  transitions  contribute  to  the  spectrum.  The  P 
branch  peak  intensity  and  its  underlying  band  head  make  this  an  interesting  case  study. 
Fig.  6  shows  the  calculated  spectrum  for  slight  changes  in  the  rotational  constants  to  B 
=  0.295  and  C  =  0.283  cm“^  for  NaH20  and  B+  =  0.298  and  C+  =  0.286  cm”^  for 
NaH2  0+  but  with  the  previous  value  for  the  A  and  A+  rotational  constants.  This  change 
reflects  about  a  2%  decrease  in  the  Na-  •  -O  bond  length  from  2.36  to  2.32  A.  These  new 
rotational  constants  result  in  a  shift  of  only  about  2  cm“^  in  the  location  of  the  peak  of 
the  rotational  contour.  However,  the  calculated  spectrum  shows  a  dramatically  different 
profile.  It  is  now  the  Q  branch  (AN  =  0)  which  has  the  highest  peak  intensity,  even  though 
it  contributes  only  13%  of  the  cross  section,  since  the  similarity  of  the  rotational  constants 
in  the  neutral  and  ion  precludes  formation  of  P  or  R  band  heads.  The  apparent  agreement 
between  these  calculated  and  measured  spectra  is  quite  good  and  may  even  suggest  that 
the  actual  rotational  constants  may  be  close  to  these  modified  values.  Comparison  of 
experimental  and  calculated  constants  for  such  clusters  can  have  large  uncertainties  due  to 
excessive  vibrational  averaging,  present  in  simple  fits  to  spectra,  but  which  would  require 
an  accurate  knowledge  of  the  global  intermolecular  potential  energy  surface  to  be  included 
in  theoretical  treatments.  In  this  context,  both  the  intermolecular  stretching  and  out-of- 
plane  bending  modes,  may  contribute  to  such  certainties  here  in  the  rotational  constants. 

Fig.  7  shows  measured  and  calculated  ZEKE  photoelectron  spectra  for  ionization  of 
the  3sai  orbital  of  the  X  ^Ai  state  of  NaNHa.  The  dominant  peak  in  the  measured  and 
calculated  spectra  again  belongs  to  P  branch  transitions  and  is  also  due  to  a  band  head. 
The  calculated  spectra  also  show  a  higher  intensity  profile  at  large  photon  energies  than 
seen  in  the  measured  spectra.  Overall  agreement  between  these  spectra,  however,  is  quite 
good. 

Comparison  of  these  spectra  with  those  of  NaH2  0  can  reveal  differences  due  to  sol¬ 
vation  effects.  For  example,  the  3sai  orbital  of  NaNHs  has  69%  s  and  24%  p  character 
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about  the  center-of-mass  compared  to  values  of  91%  s  and  7%  p  in  the  case  of  NaH20. 
This  increased  delocalization  in  NaNHs  should  lead  to  changes  in  the  photoelectron  spec¬ 
tra.  Analysis  of  the  calculated  spectrum  of  Fig.  7  shows  that  |AiV|  =  0  —  4  transitions 
contribute  14%,  45%,  23%,  10%,  and  7%  of  the  total  cross  section.  Again  the  dominance 
of  the  lAA"!  =  1  branches  is  due  to  a  surprisingly  strong  and  nonatomiclike  s-wave  com¬ 
ponent  of  the  photoionization  matrix  element.  Also,  the  close  to  20%  contribution  from 
I A  A”  I  =  3  and  4  transitions  arises  from  /-wave  (£  =  3)  components  of  the  photoelectron 
matrix  element.  These  |AA|  >  3  transitions  are  about  twice  as  large  in  NaNHs  than  in 
the  NaHaO  complex.  Furthermore,  with  their  large  Na-  •  -N  and  Na-  ■  -O  bond  lengths,  one 
would  expect  photoionization  of  these  complexes  to  occur  mainly  along  the  a  principal 
axis.  This  is  indeed  the  case  of  NaH20  where  98%  of  the  cross  section  is  associated  with 
transitions  possessing  no  change  in  angular  momentum  about  Kq,  i.e.,  AK  =  0.  However, 
with  ammonia  as  a  solvent,  only  about  60%  of  the  cross  section  is  associated  with  AKa  —  0 
transitions.  Such  differences  provide  useful  insight  into  the  structure  of  these  complexes. 

Improved  resolution  in  the  measured  ZEKE-PFI  spectra  would  provide  even  greater 
insight  into  the  structure  and  dynamics  of  these  clusters.  Several  routes  to  such  spectra 
are  being  explored.  For  example,  in  contrast  to  the  pick-up  source  currently  used,  laser 
vaporization  offers  colder  rotational  temperatures  and  would  dramatically  decrease  the 
number  of  levels  contributing  to  the  ZEKE  photoelectron  signal.  Alternatively,  multipho¬ 
ton  ZEKE-PFI  studies  can  be  performed  in  which  one  laser  selects  an  individual  eigenstate 
of  the  neutral  precursor,  which  is  subsequently  excited  to  a  high-n  Rydberg  level  and  pulse- 
field  ionized.  For  these  clusters,  the  first  step  could  involve  either  electronic  excitation  of 
the  metal  itself  or  vibrational  excitation  of  intermolecular  modes  within  the  solvent.  Such 
efforts  are  currently  underway  for  these  Na(H20)  and  Na(NH3)  complexes. 
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III.  Femtosecond  Energy-  and  Angle- Resolved  Pump-Prohe  Photoelectron 
Spectra 

1.  Background  and  Objectives 

With  femtosecond  (fs)  laser  techniques  it  is  now  possible  to  observe  nuclear  dynamics 
and  to  chart  the  path  of  chemical  reactions  in  real  time.®“^°  A  key  idea  of  this  femtosecond 
spectroscopy  is  to  use  an  ultrashort  laser  pulse,  the  pump  pulse,  to  prepare  a  coherent 
superposition  of  molecular  eigenstates  which  evolves  in  accordance  with  the  time  scales 
for  the  vibrational  (~  s)  and  rotational  (~  10"^°  s)  motions.  The  evolution  of  the 

wavepacket  is  then  probed  by  a  time-delayed,  second  ultrashort  laser  pulse  via  excitation  to 
a  final  state  which  serves  as  a  screen  for  projection  of  its  dynamics.  A  variety  of  detection 
techniques  using  different  types  of  final  states  has  been  used  to  probe  the  ultrafast  dynamics 
of  these  wavepackets.®’®’^^  These  include  the  use  of  absorption,  laser-induced  fluorescence 
(LIF),  multiphoton  ionization  (MPI),  photoelectron  spectroscopy,  nonlinear  degenerate 
four-wave  mixing,  time-resolved  mass  spectroscopy,  and  stimulated  emission  pumping. 

Of  these  techniques  the  use  of  ionization  for  detection  in  femtosecond  pump-probe 
experiments  offers  several  advantages.  The  ground  state  of  an  ion  is  stable  and  is  often  the 
best  characterized  excited  state  of  a  molecule.  Furthermore,  detection  of  ions  provides  mass 
and  kinetic-energy  resolution  in  time-resolved  photodissociation  experiments.^^  Ionization 
also  provides  photoelectrons  and  measurements  of  their  spectra  as  a  function  of  delay  time 
between  the  pump  and  probe  pulses  can  provide  additional  and  complementary  informa¬ 
tion  on  the  evolution  of  the  wavepacket  in  real  time.  Engel  and  coworkers^^  have  illustrated 
how  well  suited  photoelectron  kinetic  energy  distributions  from  pump-probe  ionization  at 
different  delay  times  are  for  mapping  wavepacket  motion  in  the  diatomic  systems  Na2  and 
Nal.  The  utility  of  this  technique  has  been  demonstrated  experimentally  on  electronically 
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excited,  states  of  sodium  dimer. More  recently,  wavepacket  evolution  on  the  A  excited  s- 
tate  of  Nal  has  been  monitored  by  detecting  both  photoions  and  photoelectrons. Charron 
and  Suzor-Weiner  have  also  presented  a  wavepacket  study  of  the  femtosecond  dynamics  of 
this  system.^® 

While  these  studies  have  clearly  served  to  illustrate  the  utility  and  promise  of  pho¬ 
toelectron  spectroscopy  for  real-time  mapping  of  wavepacket  dynamics  in  small  molecular 
systems,  they  have  generally  ignored  any  dependence  of  the  underlying  photoionization 
amplitudes  on  geometry  and  energy  in  their  studies  of  the  energy  distributions  of  the 
photoelectrons  and  have  also  not  examined  the  associated  angular  distributions  of  these 
spectra.  Neglecting  the  dependence  of  the  photoionization  matrix  element  on  geometry 
can  generally  be  a  poor  approximation  if  the  wavepacket  moves  through  a  region  of  an 
avoided  crossing  where  the  electronic  character  of  the  wavefunction  can  change  dramati¬ 
cally  or  samples  a  large  range  of  internuclear  distances.  In  fact,  the  double-minimum  state 
of  Naa  studied  by  Engel  and  coworkers^^  and  the  state  of  Nal  studied  by  Charron  and 
Suzor- Weiner are  both  formed  by  an  avoided  crossing  of  two  diabatic  states.  Further¬ 
more,  photoelectron  angular  distributions  can  be  expected  to  convey  richer  structural  and 
dynamical  information  than  is  contained  in  angle- averaged  energy  spectra. 

The  objective  of  our  effort  has  been  to  exploit  femtosecond  energy-  and  angle-resolved 
photoelectron  spectra  as  a  probe  of  wavepacket  dynamics.  We  have  developed  a  formula¬ 
tion  of  time-resolved  photoelectron  spectroscopy  that  is  well  suited  for  such  applications 
and  have  used  it  to  study  wavepacket  motion  in  the  ^  double-minimum  state  of  Na2.  In 
these  studies  we  use  geometry-  and  energy-dependent  photoionization  amplitudes  derived 
from  calculations  employing  sophisticated  descriptions  of  the  wavefunctions  for  the  double¬ 
minimum  state  and  for  the  molecular  photoelectrons  of  the  final  ionized  state.  To  illustrate 
the  utility  of  femtosecond  photoelectron  spectra  for  probing  both  molecular  vibrations  and 
rotations  in  real  time,  we  examine  spectra  for  molecules  aligned  via  a  linearly  polarized 
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pump  pulse  and  photoionized  by  a  probe  pulse  polarized  either  parallel  or  perpendicular 
to  the  pump  pulse. 

These  results  illustrate  several  key  points.  First,  a  robust  description  of  the  underly¬ 
ing  photoionization  amplitudes  greatly  enhances  the  utility  of  femtosecond  photoelectron 
spectroscopy  as  a  probe  of  wavepacket  dynamics.  In  fact,  the  use  of  such  photoionization 
amplitudes  is  essential  if  the  wavepacket  moves  through  an  avoided  crossing  or  samples 
a  large  region  of  internuclear  distances.  Second,  the  angular  distributions  of  the  photo¬ 
electrons  are  insightful  fingerprints  of  vibrational  wavepacket  dynamics.  Furthermore,  the 
strong  dependence  of  these  angular  distributions  on  the  relative  polarizations  of  the  pump 
and  probe  pulses  can  be  readily  exploited  to  monitor  molecular  alignment. 

2.  Highlights  of  Progress 
(a)  Formulation 

Fig.  8  illustrates  key  features  of  the  scheme  used  in  our  studies  of  energy-  and  angle- 
resolved  photoelectron  spectra  for  femtosecond  ionization  of  vibrational  wavepackets  in  the 

double- minimum  state  of  aligned  Na2  molecules.  A  femtosecond  polarized  laser  pulse 
of  frequency  oji  prepares  a  wavepacket  on  the  double- minimum  state  which  is  then 

ionized  by  a  time-delayed  polarized  femtosecond  laser  pulse  of  frequency  U2.  The  polar¬ 
ized  pump  pulse  produces  an  aligned  distribution  of  Na2  molecules  {Or)  since  only  those 
molecules  with  their  transition  moments  fl  parallel  or  nearly  parallel  to  the  polarization 
vector  of  the  pump  pulse  (cpump)  are  excited.  The  vibrational  and  rotational  motions  of 
this  system  with  their  respective  time  scales  of  10“^^  and  10~^°  seconds  are  monitored  via 
the  energy  and  angular  {0k,4>k)  distributions  of  photoelectrons  produced  from  ionization 
of  the  wavepacket  by  a  pulse  with  a  polarization  vector  ^ump(^P)  <^p)-  Both  tpump  and 
^probe  lie  in  the  yz-plane. 
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The  time-dependent  wavefunction  for  this  system  can  be  expanded  as 


'if  (^r,  R,  =  xo  (R.  t)  $0  (r;  R)  +  Xe  [R,  t)  (r;  R)  +  J  ^  R)  >  (2) 

where  $0)  respective  eigenfunctions  for  the  ground,  excited,  and  final 

ionized  states,  respectively,  xo,  Xe,  and  Xk  nuclear  wavepackets  on  the  individual 

potential  energy  surfaces,  r  denotes  electronic  coordinates,  and  R  is  the  internuclear  vector. 
The  interaction  between  the  molecule  and  the  laser  fields,  V{t)^  is  given  by  (see  fig.  8) 

V{t)  =  Vi{t;6R,(j)R,u!i)  +  V2{t]dR,(l)R,0p,u;2]^T) 

=  Eifi{t)sm{ujit)ep  ump  P  +  2E2f2{t  -  AT)  cos{u}2t)eprobe  '  P,  (3) 

where  Ei  and  E2  are  the  electric  field  amplitudes  of  the  pump  and  probe  pulses,  fi  (t)  and 
f2{i)  are  the  envelopes  of  these  pulses,  fj,  is  the  electric-dipole  operator,  and  AT  is  the 
time  delay  between  the  pump  and  probe  pulses.  With  ^(t)  and  V{t)  of  Eqs.  (2)  and  (3), 
projection  of  ($o(-R)|,  {^e{R)\,  and  onto  the  Schrodinger  equation  yields  equations 

of  motion  for  the  nuclear  wavepackets  xo(-R,^),  Xe(-R,^),  and  x^(-R,0 

ih^XoiR,t)  =  [Tn  +  Vo{R)]xo{R,t)  +  {^oiR)\yi{t',^RAR,^i)\^e{R))Xe{R,t)  (4) 

ih-^Xe{Rj)  =  [Tn  +  Ve{R)]XeiR,'t)  +  {^e{R)\Vl{t-,  Or,  <j)R,  U)i)\^o{R))XoiR,t) 

+  J  <?!>«, ^P,^2;Ar)|$^")(^))x,-(^,0  (5) 

and 

a  (khy 

ih^Xi(RA  = 

where  14 (i?),  14 (-R),  and  Vion{R)  are  the  potential  energy  surfaces  for  the  ground,  excited, 
and  ionic  states,  respectively,  and  m  is  the  mass  of  the  electron. 


15 


For  a  linearly  polarized  pulse,  the  interaction  matrix  element  between  the  ground  and 
excited  states  is  given  by 

{^e\Vi{t;9R,(t)R,u}i)\^o)  =  //eo  cos(6'/j)  (7) 

where  fieO  is  the  electric-dipole  matrix  element  between  the  ground  and  excited  states. 
To  obtain  the  coupling  matrix  element  between  the  excited  and  ionized  states,  we  write 
R)  as  antisymmetrized  product  of  an  ion  wavefunction  (Na^)  and  a  photoelectron 
orbital  With  an  expansion  of  <^1”^  in  its  partial  wave  components,  ,  the  coupling 

matrix  element  between  the  excited  state  $e  the  final  ionized  state  can  be  written  as 

{^^r\R)\V2{t-,eR,<t>R,dp,u2-,AT)\^e{R)) 

=  2E2  •  /2(t  -  AT)  cos(a;2t)  ^  CemYtm{k)  (8) 

im 

These  Cem  coefficients,  which  are  obtained  using  numerical  solutions  of  the  partial  wave 
components  of  the  photoelectron  orbital,  provide  the  underlying  dynamical  informa¬ 

tion  for  photoionization  of  an  oriented  Na2  molecule. 

With  a  similar  expansion  of  the  ion  nuclear  wavepacket  partial  waves  about  k, 

the  direction  of  the  photoelectron, 

==^Xklm  (^R,tjYlrn  (k'j  ,  (9) 

Im 

Eqs.  (5)  and  (6)  can  be  written  as 

ih—Xe{R,i)  =  [^AT  +  Ve]Xe{R,t)  +  VeOXo(-R,t) 

-t-  /  dkk'^'^E2f2{t  -  ^T)exp{iLO2t)Ci^{k,6R,(f)R,6p)xk£m{R,t)i'i-0) 
im 

and 

d  r  (kh)^'\ 

ikl-^Xklm{,RtY}  =  Tl<[  "t”  Vxon  2m~  Xktm.{R')^) 

+  E2f2{t  -  AT)  exp{-iuj2t)Cem{k,  9r,  (f>R,  9p)xe{R,t)-  (11) 
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where  in  Eq.  (10) 

VeO  —  Voe  =  ($e|Vi  {t;  9R,(j)R,UJi)  l$o)-  (12) 

Discretization  of  the  continuum  integration  in  Eq.  (10)  via  a  Gaussian  quadrature  of  order 
Nfc  leads  to  coupled  equations  for  Xo(-?2,^),  Xe(-R,<),  and  the  set  of  partial-wave  functions 
{Xkjtm}  for  each  quadrature  point  kj.  The  coupled  equations  are  solved  with  a  split- 
operator  technique.  In  the  present  applications  use  of  14  Gaussian  quadrature  points 
provided  converged  resxolts. 

(b)  Potential  curves  and  dipole  matrix  elements 

The  potential  energy  curves  for  the  ground  states  ofNa2  {X  ^S^)  andofNa^  (X  ^E+) 
were  obtained  at  19  internuclear  distances  between  2  and  11  A  from  extensive  Cl  calcula¬ 
tions.  Extensive  Cl  wavefunctions  were  used  for  the  ground  and  double-minimum  states 
in  the  calculation  of  the  dipole-matrix  elements  for  excitation  of  the  ground  state  and  for 
photoionization  of  the  excited  state.  Such  a  wavefunction  is  essential  to  describe  the 
state  of  Na2  which  results  from  an  avoided  crossing  of  two  adiabatic  states. The  first  of 
these  gives  rise  to  the  inner  well  and  is  a  Rydberg  state  with  about  96  %  p  character.  The 
outer  well  arises  from  the  second  diabatic  state  which  has  considerable  ionic  character  at 
large  distance  but  whose  description  changes  with  nuclear  separation.  This  behavior  with 
nuclear  separation  results  in  a  dipole-excitation  amplitude  that  is  very  Franck-Condon  in 
the  inner  well  but  changes  substantially  in  the  outer  well.  The  very  molecular  nature  of 
this  excited  state  at  very  large  internuclear  distances  can  be  seen  in  the  shown  in 
Fig.  8. 

To  obtain  the  photoionization  amplitudes  we  again  employed  the  Cl  wavefunction  for 
the  double-minimum  excited  state  but  used  a  frozen-core  Hartree-Fock  description  of  the 
wavefunction  for  the  ionized  state.  In  this  model  the  wavefunction  is  an  antisymmetrized 
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product  of  Hurtrcc-Fock  ion  orbitnls  mid  3>  photodcctron  orbitnl  thnt  is  n  solution  of  n  one 
electron  Schrodinger  equation  containing  the  Hartree-Fock  potential  of  the  Na^  ion.  We 
have  previously  developed  numerical  procedures  used  to  obtain  the  photoelectron  orbitals.^ 
The  Hartree-Fock  wavefunction  provides  a  very  adequate  description  of  the  Naj”  system 
over  the  complete  range  of  internuclear  distances  of  interest  here. 

The  changing  character  of  the  excited  state  wavefunction  with  geometry  is  expected 
to  result  in  a  significant  dependence  of  the  photoelectron  matrix  elements  and  hence  of 
the  dynamical  coefficients  Ctm  of  Eq.  (8)  on  internuclear  distance.  In  fact,  this  results  in 
an  unusual  oscillatory  behavior  of  the  Ctm  coefficients  with  geometry  around  the  barrier 
and  across  the  outer  well.  The  coefficients  for  a  photoelectron  energy  of  0.5967  eV 

and  for  the  polarization  vectors  of  the  pump  and  probe  pulses  parallel  to  the  molecular 
axis  (only  m  =  0  terms  allowed)  shoAvn  in  Fig.  8  illustrate  this  behavior  well. 

(c)  Non-Franck- Condon  behavior  and  oscillation  of  ion  signal 

A  significant  objective  of  this  work  is  to  highlight  the  importance  of  incorporating 
realistic  and  geometry-dependent  photoionization  amplitudes  in  calculations  of  femtosec¬ 
ond  photoelectron  spectra  for  wavepacket  motion.  An  excellent  illustration  of  this  arises 
from  noting  that  all  \Cern{kj^R)\'^  for  all  relevant  {kj,£,m)  essentially  vanish  in  a  narrow 
range  to  the  left  of  the  potential  barrier  on  the  excited  state  4  A).  This  suggests  that 
ionization  of  a  wavepacket  placed  in  the  inner  well  by  the  pump  laser  should  result  in  a 
very  low  photoelectron  signal  whenever  the  wavepacket  hits  this  “depleted”  region  and 
hence  in  an  oscillatory  signal  as  a  function  of  delay  time.  This  behavior  would  be  highly 
unusual. 

To  illustrate  this  feature  we  select  a  pump  photon  energy  of  3.62  eV  which  is  0.06  eV 
below  the  top  of  the  barrier.  The  lowest  vibrational  level  of  the  ground  state  is  pumped 
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with  a  FWHM  120  fs  Gaussian  pulse  centered  at  t  =  0,  forming  a  wavepacket  around  3.1  A 
directly  above  the  initial  wavepacket.  This  wavepacket  immediately  begins  an  oscillatory 
motion  with  a  period  of  347  fs,  striking  the  outer  turning  point  around  4.1  A,  where  the 
ion  signal  should  be  very  weak,  at  158  fs.  The  effect  of  the  pump  pulse  is  negligible  after 
t  -  145  fs  and  hence  the  excited  state  population  is  constant  beyond  this  time.  This  excited 
state  wavepacket  is  ionized  by  a  laser  pulse  with  a  FWHM  of  40  fs,  a  photon  energy  of  2.28 
eV,  and  polarization  parallel  to  that  of  the  pump  pulse.  The  total  ion  signal,  obtained  as 
an  integral  over  all  photoelectron  energies  and  angles,  i.e.. 


=  J  dQkJ  dR  \x^{R,tf)\'^  =  jdkk'^^JdR  \xkem{RJf)\\  (13) 


is  shown  for  several  delay  times  between  t  —  121  fs  and  1270  fs  in  Fig.  2(a).  tf  in  Eq.  (13) 
is  taken  sufficiently  long  after  the  probe  pulse  interaction  is  over.  The  ion  signal  exhibits 
an  oscillatory  behavior  with  a  period  determined  by  the  motion  of  the  wavepacket  in  the 
inner  well. 

To  provide  further  insight  into  these  minima  in  the  ion  signal,  we  show  the  form  of  the 
wavepacket  at  times  of  610,  726,  and  842  fs.  These  times  are  indicated  by  (I),  (II),  and  (III) 
in  Fig.  (9).  At  610  fs  the  wavepacket  is  delocalized  and  the  ion  signal  results  from  ionization 
over  a  wide  region  of  this  inner  well  ((I)  in  Fig.  9).  After  passing  through  the  intermediate 
stage  at  726  fs  (II),  there  is  a  significant  accumulation  of  the  wavepacket  around  4  A,  the 
righthand  turning  point,  at  842  fs  (III).  Depletion  of  the  dynamical  coefficients  Ctm  results 
in  a  dramatic  decrease  in  the  ion  signal  for  ionization  at  this  time.  A  robust  description 
of  the  photoionization  dynamics  is  clearly  essential  for  capturing  this  behavior. 


(d)  Photoelectron  energy  distributions 


We  next  look  at  the  photoelectron  energy  distributions  for  ionization  of  a  wavepacket 
that  is  formed  by  a  pump  pulse  with  an  energy  almost  equal  to  the  top  of  the  poten¬ 
tial  barrier  in  the  excited  state  {huji  =  3.68  eV).  Fig.  10  shows  the  behavior  of  such  a 
wavepacket  formed  by  a  pulse  with  a  FWHk/I  of  120  fs  and  centered  at  t  0.  It  takes  the 
wavepacket  about  100  fs  to  reach  the  barrier  at  4.7  A  where  it  splits  into  two  components, 
a  lower-energy  component  that  is  reflected  into  the  inner  well  and  reaches  its  inner  turning 
point  at  <  =  400  fs  and  800  fs  and  a  higher-energy  component  in  the  outer  well  that  reaches 
its  outer  turning  point  at  t  =  600  fs  and  rejoins  the  inner  well  component  at  t  =  1000 
fs.  This  wavepacket  is  photoionized  by  a  probe  pulse  with  energy  hu2  of  2.28  eV  and  a 
FWHM  of  40  fs.  The  photoelectron  signal  in  an  energy  interval  dek  is  given  by 

Pion{ek)dek  =  kdtk  f  dQk  f  dR  =  kdtkY,  /  dR  \xkim{R.i  f)?  ■  (14) 

J  J  im  d 

Fig.  11  shows  the  resulting  photoelectron  energy  distributions  as  a  function  of  delay 
time  AT  between  the  pump  and  probe  pulses  and  with  the  polarization  vectors  of  the  pump 
and  probe  parallel.  Characterization  of  the  peaks  in  this  spectrum  is  straightforward. 
Fig.  11  shows  that  the  wavepacket  moves  slowly  through  the  barrier  region  at  around 
AT  =  200  and  1000  fs.  The  large  values  of  the  dynamical  coefficients  sandwiching 

this  barrier  region  between  4  and  5  A  (see  Fig.  8)  result  in  strong  photoelectron  peaks 
around  0.6  eV  for  these  delay  times  of  200  and  1000  fs.  A  strong  peak  at  very  low 
photoelectron  energy  (~  0.0002  eV)  is  seen  around  AT  =  600  fs  when  the  wavepacket 
hits  the  outer  turning  point  near  R  =  S.b  k  where  only  photoelectrons  with  low  kinetic 
energy  can  be  ejected  since  almost  all  the  energy  goes  into  potential  energy  of  NaJ".  While 
the  signal  is  less  intense  for  ionization  of  this  wavepacket  by  a  probe  pulse  polarized 
perpendicular  to  the  pump  pulse,  the  photoelectron  spectrum  is  very  similar  to  that  of 

Fig.  11. 
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(e)  Photoelectron  angular  distributions 

Angular  distributions  are  well-known  to  provide  valuable  insight  into  the  underlying 
dynamics  of  conventional  photoelectron  spectra.  In  time-resolved  pump-probe  photoelec¬ 
tron  spectroscopy,  these  angular  distributions  can  also  be  expected  to  be  significant  dy¬ 
namical  windows  on  wavepacket  behavior  and  on  the  evolving  electronic  structure.  We 
now  report  the  first  results  of  such  studies.  These  time-domain  photoelectron  angular 
distributions  are  given  by 

A[6k')  ^k)d^k  —  j  dR  ^'^)^im{^kj '^)\  kdef;^  (lb) 

d  tm 

where  6k  is  measured  relative  to  the  polarization  of  the  probe  pulse.  Fig.  12  shows  these 
angular  distributions  for  the  probe  pulse  (2.28  eV)  polarized  parallel  and  perpendicular 
to  the  pump  pulse  (3.68  eV)  and  for  a  photoelectron  energy  of  0.5967  eV.  In  contrast 
to  the  similarity  of  the  photoelectron  energy  distributions,  the  angular  distributions  are 
strikingly  different  for  the  “parallel”  and  “perpendicular”  cases.  For  the  parallel  case 
the  angular  distributions  are  clearly  of  the  3.^2  type,  as  shown  in  the  inset  of  Fig.  12, 
while  they  show  dyz  behavior  for  the  perpendicular  case.  This  behavior  is  consistent 
with  symmetry  considerations  for  the  dipole  interaction.  Fig.  10  shows  that  the  weaker 
photoelectron  signal  for  ionization  at  605  fs  comes  from  the  inner  well  region  where  the 
photoionization  amplitude  is  low  while  the  stronger  signal  at  AT  of  968  fs  comes  from 
the  barrier  region  as  the  wavepacket  returns  from  the  outer  basin.  In  this  barrer  region 
the  wavepacket  moves  slowly  and  the  ionization  amplitudes  are  laxge.  The  corresponding 
partial  wave  components  of  the  ion  wavepacket  Xkjtm{R,t)  provide  very  graphic  images  of 
the  photoionization  dynamics  of  these  wavepackets.  Some  examples  are  shown  in  Fig.  13. 

The  strong  dependence  of  the  photoelectron  angular  distributions  on  the  angle  between 
the  molecular  axis  and  the  polarization  of  the  probe  pulse  makes  these  distributions  a 
valuable  real-time  probe  of  molecular  rotation.  In  the  present  example,  consider  a  molecule 
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rotating  in  the  yz-plane  and.  oriented  at  an  angle  relative  to  the  z  axis  at  a  probe  time 
AT.  For  ionization  by  a  probe  pulse  polarized  along  the  z-axis,  the  angular  distributions 
in  the  direction  dk  should  behave  like 

A[Sk\^R{^T))  =  cos^ dn  X  -  distribution)  +  sin^  Or  x  {dy^  -  distribution) 

+  cos  9r  sindi?  X  (cross  term).  (16) 

The  photoelectron  signal  hence  depends  strongly  on  the  position  of  the  detector.  For 
example,  for  the  detector  along  the  z-axis,  there  should  be  a  significant  photoelectron 
signal  (due  to  the  d^2  distribution)  only  for  those  delay  times  when  the  molecules  are  lined 
up  along  the  z-axis.  Two  such  peaks  should  be  recorded  in  a  single  rotational  period. 
On  the  other  hand,  for  those  delay  times  when  the  molecule  lies  along  45°,  a  significant 
signal  should  be  detected  due  to  the  dj,^-type  distributions.  Four  such  large  signals  should 
be  observed  in  a  single  period.  With  the  detector  along  the  y-axis,  only  the  dj^^-type 
distributions  would  result  in  a  significant  photoelectron  signal.  This  peak  signal  should 
occur  for  0r{AT)  w  45°.  No  significant  signal  is  expected  from  the  d^2-type  distributions. 
This  strong  dependence  of  the  photoelectron  signal  on  delay  time  can  be  used  for  real-time 
mapping  of  molecular  rotation  and  can  be  a  valuable  supplement  or  alternative  to  LIF 
detection. 

(f)  Outlook 

This  study  represents  a  significant  advance  towards  exploiting  the  full  potential  of 
femtosecond  photoelectron  spectroscopy  as  a  probe  of  wavepacket  dynamics.  The  fem¬ 
tosecond  energy-  and  angle-resolved  photoelectron  spectra  of  these  studies  constitute  very 
compelling  evidence  of  the  unprecedented  insight  that  such  spectra  can  provide  into  both 
wavepacket  motion  and  the  evolution  of  the  associated  electronic  structure.  These  time-. 
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energy-,  and.  angle-resolved  pliotoelectron  spectra  should  be  particularly  useful  for  studies 
of  wavepackets  in  the  important  cases  of  avoided  crossings  and  large-amplitude  motion. 
Because  these  studies  must  rely  on  a  robust  description  of  the  underlying  photoionization 
amplitudes,  applications  to  small  polyatomic  systems  should  be  most  rewarding.  Such  s- 
tudies  of  femtosecond  energy-  and  angle-resolved  photoelectron  spectra  in  small  polyatomic 
systems  are  underway. 
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Captions 

Fig.  1  Ion  rotational  distributions  for  (2+1')  REMPI  of  H2O  via  the  4i3  rotational  level  of 
the  C  ^B\  (Spai)  state.  Only  type  b  transitions  {5Ka  and  AKc  =  odd  or  AKa  and 
AKc  =  even)  are  allowed.  The  selection  rule  is  AKa  +  i  =  odd.  The  photoelectron 
energy  is  about  0.85  eV. 

Fig.  2  Partial  wave  components  of  the  photoionization  matrix  element  for  the  Spaj  orbital 
(«  14%  s,  84  %  p)  of  the  C  ^Bi  state  of  H2O. 

Fig.  3  The  (a)  measured  and  (b)  calculated  PFI-ZEKE  photoelectron  spectra  for  (1+1') 
REMPI  via  the  R(7)  line  of  the  A  ^IIiu  X  band.  The  e/f  labels  indicate  the 
parity  of  the  ion  rotational  levels. 

Fig.  4  The  (a)  measured  and  (b)  calculated  PFI-ZEKE  photoelectron  spectra  for  (1+1') 
REMPI  via  the  R(6)  line  of  the  A  ^n2u  band.  The  e/f  labels  indicate 

the  parity  of  the  ion  rotational  levels.  The  calculated  spectrum  is  convoluted  with 
a  Gaussian  detection  function  with  an  FWHM  of  1.75  cm“^  and  assumes  a  kinetic 
energy  of  50  meV. 

Fig.  5  (a)  Computed  ZEKE  photoelectron  spectra  for  the  ground  state  of  Na(H20)  with 
calculated  molecular  geometries  (Na-  •  -0  and  Na"^  •  •  -0  bond  lengths  of  2.355  A  and 
2.230  A,  respectively);  (b)  measured  ZEKE  spectra  of  Geoffrey  Blake  and  coworkers. 

Fig.  6  (a)  Computed  ZEKE  photoelectron  spectra  for  the  ground  state  of  Na(H20)  with 
calculated  molecular  geometries  (Na-  •  -0  and  Na"*"  •  •  -0  bond  lengths  of  2.32  A  and 
2.270  A,  respectively);  (b)  measured  ZEKE  spectra  of  Geoffrey  Blake  and  coworkers. 

Fig.  7  (a)  calculated  and  (b)  measured  ZEKE-PFI  photoelectron  spectra  for  NaNHs  at  a  ro¬ 
tational  temperature  of  100  K.  The  calculated  spectrum  assumed  rotational  constants 
calculated  at  equilibrium  geometry. 


26 


pjg  8  Potenti&l  surfaces  of  Na2-  The  transition  matrix  elements  are  superposed.  The  pho 

-  2  2 
toionization  amplitude  shown  is  defined  as  Ckjim  —  \kjy/WjCimik,  4>r-  ^p)  |  • 

Wj  is  the  quadrature  weight  at  kj.  The  scales  for  /i^e  ICfmP  3-re  to  the  right. 


Fig.  9  (a)  Oscillatory  behavior  of  the  total  ion  signal  from  a  wavepacket  confined  in  the  inner 
potential  basin,  (b)  The  spatial  distribution  of  the  wavepackets  associated  with  the 
total  ion  signals  at  the  points  (I),  (II),  and  (III)  marked  in  (a).  The  minima  of  the 
ion  signal  are  formed  by  the  packet  that  is  localized  at  the  right-hand  turning  point, 
where  the  ionization  amphtudes  are  depleted. 

Fig.  10  Spatio-temporal  distribution  of  the  wavepacket  pumped  to  the  top  of  the  barrier  on 
the  excited  state.  The  zone  between  the  two  dotted  lines  indicates  the  potential 
barrier  region  where  the  dipole-ionization  amplitudes  are  large  and  the  wavepacket 
moves  slowly. 

Fig.  11  Photoelectron  signal  versus  kinetic  energy  and  delay-time  for  a  packet  with  energy 
almost  at  the  top  of  the  potential  barrier  and  for  parallel  polarizations  of  the  pump 
and  probe  pulses. 

Fig.  12  (a)  Photoelectron  angular  distributions  as  a  function  of  the  delay  time  for  parallel 
polarizations  of  the  pump  and  probe  pulses.  The  inset  shows  a  polar-coordinate  {6^) 
description  at  delay  times  of  AT  =  605  and  968  fs.  (b)  As  in  (a)  but  for  the  probe 
pulse  polarized  perpendicular  to  the  pump  pulse. 
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